Alignment-Dependent Ionization of N2, 02, and CO2 in Intense Laser Fields 
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The ionization probability of N2, 02, and CO2 in intense laser fields is studied theoretically as 
a function of the alignment angle by solving the time-dependent Schrodinger equation numerically 
assuming only the single- active-electron approximation. The results are compared to recent experi- 
mental data [D. Pavicic et al., Phys. Rev. Lett. 98, 243001 (2007)] and good agreement is found for 
N2 and 02. For CO2 a possible explanation is provided for the failure of simplified single-active- 
electron models to reproduce the experimentally observed narrow ionization distribution. It is based 
on a field-induced coherent core-trapping effect. 
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Time-resolved imaging of the dynamics of nuclei and 
electrons on a femtosecond or even sub-femtosecond time 
scale is a prerequisite for the experimental investigation 
of the formation and breaking of chemical bonds. Ul- 
trashort laser pulses have recently been demonstrated 
to allow monitoring of nuclear motion with subfemtosec- 
ond and sub- Angstrom resolution in real time [l|43|. It 
was also experimentally demonstrated that the high- 
harmonic radiation or electrons emitted in an intense 
laser pulse may in principle reveal electronic structure 
[J, y and thus have the potential for time- resolved imag- 
ing of changes of the electronic structure in, e. g., a chem- 
ical reaction. To reach this goal it is, however, important 
to understand the relation between electronic structure 
and strong- field response in molecules. 

In [61] the alignment dependence of the strong-field ion- 
ization was measured for N2, 02, and CO2. The found 
relatively strong angular dependence of the ionization on 
the incident electric field is interpreted in terms of molec- 
ular orbital theory: the measured ionization is supposed 
to map the structure of the highest occupied molecular 
orbital (HOMO). While the angle- resolved ionization of 
N2 and 02 reproduces the symmetry of the HOMO and 
was reasonably explained using a simplified strong-field 
model like molecular Ammosov-Delone-Krainov (MO- 
ADK) theory [7|, this was not the case for CO2. 

In order to find out why simple strong-field models 
like MO-ADK or the molecular strong-field approxima- 
tion (MO-SFA) [8| fail and are unable to properly re- 
produce the CO2 results, a recently introduced approach 
(SAE-TDSE) to solve the time-dependent Schrodinger 
equation (TDSE) of a moleculewithin the single-active- 
electron approximation (SAE) [9| is adopted in this Let- 
ter to calculate the alignment dependence of ionization 
of N2, 02, and CO2. It avoids most of the assumptions 
made in the derivations of MO-ADK (a pure tunneling 
theory) or MO-SFA (that ignores the Coulomb interac- 



tion between the emitted electron and the remaining ion) 
which both neglect excited states. The reliability of the 
SAE-TDSE has been demonstrated, at least for H2, in 
i9i] by a comparison to exact calculations. In this Letter 
we show that the failure of siniple models to explain the 
experimental CO2 results in ;6] appears to be a laser- 
induced dynamical coherent coupling of the HOMO and 
the core, a phenomenon we name coherent core trapping. 

The SAE-TDSE approach [9] and its present exten- 
sion to larger molecules and arbitrary orientation is only 
briefly described here. The basic assumption is that ef- 
fectively only one single active electron interacts with the 
external laser field and the corresponding time-dependent 
Schrodinger equation is then solved numerically. The 
non-relativistic time-dependent Hamiltonian is treated 
senn-classically and separated in an unperturbed time- 
independent molecular part -ffmoi and a time-dependent 
laser-electron interaction Hint{t). A fixed molecular ge- 
ometry is used to calculate the field-free molecular elec- 
tronic structure which is presently kept invariant dur- 
ing the laser pulse. Experimental values of the equilib- 
rium distances have been adopted in the calculations. 
The time-dependent wavefunctions are expressed as a lin- 
ear combination of the solutions of the time-independent 
field- free Hamiltonian -ffmoi- The latter are Kohn-Sham 
orbitals obtained from density-functional theory using 
the exchange correlation potential LB94 [10|, since it 
yields the correct asymptotic behavior. The radial parts 
of the orbitals are expanded in a _B-spline basis con- 
sisting of one set that is defined within one large cen- 
tral sphere (with radius r'^g_-^) and additional sets de- 
fined within (non-overlapping) atom-centered spheres. 
While the latter sets assure a proper molecular descrip- 
tion of the atomic cores, the former improves the de- 
scription of molecular bonds and leads to a convenient 
"box" discretization of the electronic continuum. The 
angular parts are in both cases expressed in terms of 



spherical harmonics and molecular symmetry is fully ac- 
counted for. Converged results were obtained for radii 
r^i,^ 200.35 a.u., 160 a.u., and 200.35 a.u. and the values 
of ^max 10, 12, and 14 for N2, 02, and CO2, respectively. 

In the original SAE scheme (method A) all but a sin- 
gle active electron are frozen in the time propagation. 
The passive core orbitals are therefore removed from 
the active space in the TDSE propagation, i. e. the ac- 
tive electron is constrained to remain orthogonal to the 
field-free core. However, all electrons are driven by the 
field, and due to the unitary evolution under the time- 
dependent Hamiltonian all orbitals remain orthogonal 
at all times without the need of external orthogonality 
constraints. Thus in method B all orbitals are prop- 
agated within the complete functional space, including 
the one occupied by the field-free core. Due to the men- 
tioned orthogonality, this does not violate the Pauli prin- 
ciple and insures in fact full gauge invariance of the re- 
sults, as will be explained elsewhere. Mehod B may be 
dubbed independent-active-electron (lAE) approach and 
is in fact implicit in many model-potential calculations 
which, however, neglect orthogonality altogether. 

The inclusion of symmetries up to A = 6 where 
A = 0, 1, 2, . . . refers to cr, tt, 5, . . . symmetries was suf- 
ficient to reach convergence. The total number of states 
used for time-propagation are limited by an energy cut- 
off parameter which has been chosen between 10 a.u. and 
20 a.u. in order to yield converged results. This resulted 
in about 17,000, 22,000, and 22,000 states in the calcu- 
lations for N2, 02, and CO2, respectively. 

The time profile of the applied linearly polarized laser 
pulse is represented by a cos^ vector potential and has 
in all cases shown a wavelength of 800 nm and a dura- 
tion of 40 fs. While only results for the experimental 
peak intensities are shown, a number of calculations for 
lower peak intensities were performed and no qualitative 
changes were found for the total ion yields. Thus focal- 
volume averaging should not influence the main findings 
of this work, except one case as is discussed below. At the 
end of the time-propagation calculations (performed in 
length gauge) the ionization yield is obtained from a sum 
over the populations of the discretized continuum states. 
In the case of H2 it was demonstrated by a comparison 
to full two-electron calculations that for ionization yields 
not exceeding about 10 % the SAE-TDSE results should 
be multiplied by a factor of 2 in order to account for 
the two equivalent electrons @- Correspondingly, the 
ionization yields obtained for a given initial orbital were 
multiplied with the occupation number and orbital de- 
generacy was explicitly considered. The total ionization 
is obtained by an incoherent sum of these weighted or- 
bital yields. For the experimental parameters of [6| used 
in this work this sum remains well below 10 %, except 
for 02 where the maxima almost reach 10%. Since for 
02 the contributions from orbitals other than the HOMO 
are found to be negligible, this does, however, not change 



180 



180 




180 



15 30 45 60 75 90 

Angle (degrees) 



FIG. 1: (Color online) Left panels: Normalized ionization 
yields as a function of the angle between the electric field 
and the molecular axis: experimental data (black, [y]) and 
theoretical results (method A) for the HOMO (red dots) and 
for the total yield (green dashed line). For CO2 also the result 
of method B is shown. Right Panels: Ionization yields of 
the HOMO and some lower lying orbitals (HOMO-X) as well 
as the excitation yield (E) from the HOMO. Laser peak 
intensities: N2 (1.5 x 10^* W/cm^), 02 (1.3 x lO" W/cm^) 
and CO2 (1.1 X 10" W/cm^). 



any of the conclusions of the present work. 

Figure [T] shows a comparison of the orientation de- 
pendence of the ionization yield of N2 obtained in this 
work (method A) with the experimental data in [6|. Al- 
ready the pure HOMO (3crg) contribution agrees quite 
well with experiment, and this agreement is slightly im- 
proved, if the contributions of all orbitals are included. 
As the angle between molecular axis and field vector ap- 
proaches 90° , there is an increasing contribution from the 
HOMO-1 (ItTu). Its inclusion changes the ratio of par- 
allel to perpendicular ionization from 5.7 to 4.4 which 
then is in reasonable agreement to the experimentally 



found 3.3±0.4 [6[. This confirms the importance of the 
HOMO-1 especially for perpendicular orientation, also 
found recently for high harmonic generation in 11| . The 
contribution of the HOMO-2 (2cr„) is on the other hand 
practically negligible, though at 0° it is larger than the 
one of the HOMO-1. 

For 02 the agreement between SAE-TDSE and exper- 
iment is again good, as is seen from Fig. [TJ The small 
shift of the maxima towards the horizontal axis may be 
explained by the influence of vibrational motion, since 02 
should preferentially ionize at smaller internuclear sep- 
arations than at the equilibrium one 12|. The larger 



asymmetry between parallel and perpendicular orienta- 
tion found in the SAE-TDSE is a consequence of the fact 
that both electrons in the degenerate tt^ orbitals (tt^ , 7rp 
contribute equally in parallel orientation, but effectively 
only the tt^ electron contributes in x direction. The sym- 
metry observed in the experiment is thus somewhat sur- 
prising and it may be worthwhile to check the experi- 
mental deconvolution procedure in this respect. In fact, 
an even better agreement is found for all three consid- 
ered molecules N2, O2, and CO2, if the contributions 
from the tt^ orbitals are omitted. Compared with N2, 
the alignment dependence of the ionization yield from 
the 02 HOMO with its -Kg symmetry differs clearly, in 
accordance with expectations based on the correspond- 
ing (momentum) densities (see, e.g.. Figs. 2 and 3 in 
[13[). The contributions from the HOMO-1 (3crg) and 
HOMO-2 (l7r„) are found to be practically negligible, but 
their orbital structure is again imprinted in their ioniza- 
tion yields. For parallel alignment the Sug orbital of 02 
is, e.g., much more likely to ionize than for perpendic- 
ular alignment and it shows a very similar behavior as 
the HOMO Zag of N2. The same analogy applies for 
the l7r„ orbitals of 02 and N2. The found correspon- 
dence between orbital structure and ionization behav- 
ior for HOMO and HOMO-1 of N2 contradicts resuhs of 
a recent all-electron time-dependent density-functional- 
theory calculation 



IJ]. Noteworthy, the present SAE- 



TDSE results for both N2 and 02 agree better to exper- 
iment than the results in 14| . 

A similar correlation with the orbital shape is found 
for the orientation dependence of the HOMO and the 
HOMO-X (X=l, 2, 3) of CO2. Although the HOMO 
of 02 and CO2 are both lUg orbitals, the one of CO2 
is more stretched along the nuclear axis which leads to 
a more asymmetric orientation dependence (peaking at 
about 40°). In that sense the present SAE-TDSE results 
confirm the conjecture that the alignment dependence of 
the strong-field ionization reflects the symmetry of the 
orbital. A prerequisite for this correspondence is the ab- 
sence of intermediate resonances. The excitation yields 
to all electronically excited bound states are also shown 
in Fig. [1] For all three molecules they are much smaller 
than the ionization yield and, more importantly, follow 
in shape quite well the angular dependence of the HOMO 
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FIG. 2: (Color online) Normalized ionization yield from the 
HOMO of CO2 (method B) for different laser peak intensities 
between 3 x 10^"^ W/cm^ (as in the experiment 18]) and 1.1 x 
10" W/cm^ (experiment in @|). 



ionization yields. Therefore, though excitation may not 
necessarily be negligible in absolute magnitude for the 
three considered molecules and laser parameters, it does 
not modify the orientation dependence of the ion yields. 

Nevertheless, the agreement of the CO2 result with 
experiment is not good. The latter found very sharp 
maxima, the first peaking at about 46°, and extremely 
small ionization in parallel and perpendicular orienta- 
tion. Since the contribution from the HOMO-2 (3cr„) 
is non-negligible for the parallel orientation according to 
SAE-TDSE (in agreement with recent experimental evi- 
dence [15]), the agreement gets even worse for the total 
ionization. Note, however, that the present results agree 
reasonably well with the measured angular distribution 
of fragments of CO2 in jl6!|. Furthermore, already in 
6] the experimental CO2 distribution was found to be 
disagreeing to simplified SAE models like MO-ADK and 
MO-SFA. Since the present SAE-TDSE is a stringent test 
of the SAE beyond such simplified models, it might be 
concluded that the deviation to experiment is a break- 
down of the SAE itself, if it is not experimentally caused 
as is pro posed in a recent compilation of theoretical re- 
sults [l7|. The most likely explanation for the observed 
sharp angular distribution appears to be an intermediate 
resonance whose position is either badly or (in the case of 
doubly-excited states) not at all predicted within SAE. 
Clearly, an experimental wavelength (or intensity) scan 
should reveal such a resonance and would thus be very 
desirable. In fact, there exists one further experimen- 
tal investigation of the orientation dependent ionization 
of CO2. In [18| data are presented for an about 3 to 
4 times lower peak intensity than in [(|. The distribu- 
tion is much broader and thus in better agreement with 
simplified models. 

So far, all discussed results were obtained with method 
A. While there is no qualitative change when applying 
method B to N2 and O2, this is not the case for CO2 
(see Fig.[T|). In this case the SAE-TDSE results are more 
similar to the experiment, since the peak shifts towards 



a larger angle and becomes narrower. Figure [2] shows the 
results obtained with method B for different laser intensi- 
ties and reveals an even more interesting effect. The peak 
widths become pronouncedly narrower with increasing 
intensity. This is in accordance with the experimental re- 
sults of [6| and 18|. In fact, there is a pronounced change 
occurring at a peak intensity around 8 x 10^'^ W/cm^ and 
thus in between the two experimentally investigated in- 
tensities. A more detailed analysis reveals that this effect 
is due to a field-induced coupling between the field-free 
HOMO (iTTg) and the HOMO-1 (l7r„). In fact, at an in- 
tensity of about 8 X 10^^ W/cm^ the different ac Stark 
shifts of the two involved orbitals allow strong Rabi os- 
cillations to occur and the resulting coherent trapping of 
the electron reduces the ionization probability. 

Due to symmetry, the Itt^ HOMO and the ItTu 
HOMO-1 are coupled by a field parallel to the molec- 
ular axis. Therefore, the effective intensity responsible 
for a coupling of these states is the component parallel 
to the molecular axis. For a given intensity the transition 
may be resonant for the parallel orientation, but is not 
yet resonant for other orientations. A further increase of 
the intensity allows then a larger range of orientations 
to become resonant. As a consequence, the orientational 
distribution becomes narrower, as the ionization proba- 
bility decreases for an increasing range of angles close to 
0° (or 180°). It was confirmed that the shift of the max- 
imum of the yield agrees with the corresponding shift of 
the projection of the laser field onto the molecular axis. 

Note, there are two effects of the strong field that lead 
to coherent core trapping. First, the Stark shift brings 
the orbitals into resonance. (Note, however, that for N2 
the ac Stark shift suppresses the one-photon coupling be- 
tween HOMO and HOMO-1 that in the field- free situa- 
tion would be almost resonant for 800 nm photons.) Sec- 
ond, the field-induced quiver motion temporarily depletes 
the field-free core orbital. Clearly, the present theoreti- 
cal results alone would not be conclusive. Although the 
adopted approximation of an independent response of the 
electrons should be reasonable for intense fields (cf. [y]), 
the results cannot be expected to be quantitatively cor- 
rect. In fact, due to the intensity dependence also focal- 
volume effects are relevant. However, together with the 
experimental findings the present results strongly sug- 
gest that coherent core trapping is responsible for the 
unexpected experimental results observed for CO2. The 
field-free excited states (and their population) are on the 
other hand not found to lead to substantial effects for 
the considered examples, in contrast to the interpreta- 
tion given in [19| . 

In summary, we have reproduced satisfactorily exist- 
ing experimental data for the angle-resolved ionization of 
the diatomic molecules N2 and 02 in intense short laser 
pulses. This shows the potential of the present SAE- 
TDSE approach that had so far only been tested for H2. 
These results support the principle idea of orbital map- 



ping, but we also observe important contributions from 
the HOMO-X that may distort the tomographical pic- 
ture of the HOMO as in N2. On the other hand, CO2 
is an example for an intensity-dependent breakdown of 
the orbital tomography that appears to be caused by a 
field-induced coupling of the HOMO with the core. This 
may explain why two experiments performed with differ- 
ent laser intensities observed different angular ionization 
patterns. 
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